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One designs a Stratonovich noise feedback controller with support in an arbitrary open subset Oo of O which exponentially stabilizes 
in probability, that is with probability one, the Oseen-Stokes systems in a domain O c M. d , d = 2,3. This completes the stabilization 
results from the author's work [6] which is concerned with design of an Ito noise stabilizing controller. 
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1. Introduction and statement of the problem 

This work is concerned with internal stabilization via Stra- 
tonovich noise feedback controller of Oseen-Stokes system 



dX 

vAX + (/ • V)X + (X ■ V)g = Vp in (0, oo) x O, 

at 

V • X = in (0, oo) x O, 

X(0, g) = £ € O, X = on (0, oo) x 80. 



(1) 



o 



X 



Here, O is an open and bounded subset of Mr, d = 2,3, with 
smooth boundary 80 and /, g e C 2 (0; Mr) are given functions. 
In the special case g = 0, system (1) describes the dynamic of 
a fluid Stokes flow with partial inclusion of convection acce- 
leration (/ • V)X (X is the velocity field). The same equation 
describes the disturbance flow induced by a moving body with 
velocity / through the fluid. Should we mention also that in the 
special case f = g = X e , where X e is the equilibrium (steady- 
state) solution of the Navier-Stokes equation 



ax _ 

dt 

V X 



vAX + (X ■ V)X = Vp + f e , 
= 0, X\ d0 = 0, 



(2) 



and f e e C(0; M. d ), system (1) is the linearization of (2) around 
X e . In this way, the stabilization of (1) can be interpreted as the 
first order stabilization procedure of steady-state Navier-Stokes 
flows. 

Our aim here is to design a stochastic feedback controller 
of the form 

M 

u = \ 0o J] R k {X) ofi k , R k e L{{L\0)) d ), (3) 

k=l 

which stabilizes in probability system (1) and has the support 
in an open subdomain Oo of O. 



Here, {fik}^ is a system of mutually independent Brownian 
motions on a probability space {Q, T , P) with filtration {T t }t>o 
while the corresponding closed loop system 



dX - vAXdt + (/ • V)Xdt + (X ■ V)gdt 



M 



l 0o 2^MX)odf3 k + Vpdt 

k=\ 



(4) 



X(0) 



is taken in Stratonovich sense (see, e.g., [1]) and this is the sig- 
nificance of the symbol Rk(X)o0 k in the expression of the noise 
controller (3). In the following, the terminology Stratonovich 
feedback controller or Ito feedback controller refer to the sense 
in which the corresponding stochastic equation (4) is consi- 
dered: in the Stratonovich sense or, respectively, the Ito sense. 
We have denoted by \q b the characteristic function of the open 
set Oo c O. 

In [5], [6], [7], [8], the author has designed similar stabili- 
zable Ito noise controllers for equation (1) and related Navier- 
Stokes equations. However, it should be said that, with respect 
to Ito noise controllers, the Stratonovich feedback controller (3) 
has the advantage to be stable with respect to smooth changes 
F k ' of the noise ji k and this fact is crucial not only from the con- 
ceptual point of view, but also for numerical simulations and 
practical implementation into system (1) of the random stabi- 
lizable feedback controller 



« £ (0 = l Ol Yw))^), 
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k=l 

where is a smooth approximation of /3 s k . If instead (3) we 
take u to be an Ito stabilizable feedback controller, then the 
corresponding Ito stochastic closed loop equation is convergent 
for s — > to a Stratonovich equation of the form (4) which 
might be unstable because the stabilization effect of the noise 
controller is given by the Ito's formula which is valid for Ito 
stochastic equation only. Thus, for numerical implementations 
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of a stabilizanble noise feedback controller u it is essential that 
it is of Stratonovich type. It should be said, however, that the 
option for an Ito stabilizable noise controller, like in [5]-[8], or 
a Stratonovich one, as in this work, is function of the specific 
techniques used to insert such a noise controller into system (1); 
by direct simulation or by numerical approximating procedure. 

As regards the literature on stabilization of linear differen- 
tial systems by Stratonovich noise, the pioneering works [2], [3] 
should be primarily cited. For linear PDEs, this procedure was 
developed in [9], [10] which are related to this work. For ge- 
neral results on internal stabilization of Navier-Stokes systems 
with deterministic feedback controllers, we refer to [4]. (See 
also [11] for noise stabilization effect into a different nonlinear 
PDE context.) 

2. The noise stabilizing feedback controller 

Consider the standard space of free divergence vectors 
H = [y e (L 2 (0)) d ; V • y = in O, y ■ n = on 80} and denote 
by j?lo : D(Jlo) c H -> H the realization of the Oseen-Stokes 
operator in this space, that is, 

& y = P(-vAy + (f-V)y + (yV)g), y e D(^\ (5) 

where D(&o) =Hn (H 2 (0)) d n (H l a (0)) d . Here, P is the Leray 
projector on H and H 2 (0), H^(0) are standard Sobolev spaces 
on O. In the following, it will be more convenient to repre- 
sent equation (1) in the complex Hilbert space <H = H + iH 
by extending yfo to 3i : D(3\) •z'H — > "H via standard proce- 
dure, J[(x + iy) = J[qx + i&oy. The operator has a countable 
set of eigenvalues {Aj}J =l (eventually complex) with the corres- 
ponding eigenvectors tpj. Denote by !K the adjoint of J{ with 
eigenvalues Aj and eigenvectors <p*. Each eigenvalue is repeated 
in the following according to its algebraic multiplicity mj. Nor- 
malizing the system {<fij}J =l , we see that 

l^ylw = *j " ((/ ' + (jpj ■ V)g, <pj) , V/, 

and since, by the Fredholm-Riesz theory, \Aj\ — > +oo as j — » oo, 
we infer that Re Aj — » +oo as j — » oo. We denote by the 
minimal number of eigenvalues Aj for which 

Re Aj > for ; > N, Ai + A 2 + ■ ■ ■ + A N > 0. (6) 

(In the above sequence, each Aj is taken together its conjugate 
Aj and, clearly, there is such a natural number N.) 

Set A"„=lin spanf^}^ and denote by X s the algebraic com- 
plement of X u in X. It is well known that X u and X s are both 
invariant for J{ and, if we set 

■ftu = &\x u , &s = &\x s , 

we have for their spectra cr(Jl u ) = {Aj}^, cr(Jl s ) = {Aj}°° =N+1 
and, since -ft. s is the generator of an analytic C-semigroup 
e -^,t m x s , we have 

\\e- ms, \\im <Cexp(-ReA w+ iO, t>0, (7) 
(see, e.g., [5], p. 14). In the following, we shall assume that 



(i) All the eigenvalues {Aj} . =l are semisimple. 

This means that the algebraic multiplicity of each Aj, j - l,...,N, 
coincides with its geometric multiplicity or, in other words, the 
finite-dimensional operator (matrix) JK U is diagonalizable. As 
we will see later on, this assumption is not essentially neces- 
sary but it simplifies however the construction of the stabilizing 
controller because it reduces the unstable part of the system to 
a diagonal finite-dimensional differential system. In particular, 
it follows by (i) that we can choose the dual systems {<pj} and 
[ip*\ in such a way that 

(<p i ,<p)) = 6 ij , /•-./- 1 V. (8) 

(Here, and everywhere in the following, (•, •) stands for the 
scalar product in 'H and H. By | • \<n and | ■ |# we denote the 
corresponding norms.) 

We note that the uncontrolled Oseen-Stokes system (1) can 
be rewritten in the space as 

^ + = 0, t > 0, X(0) = x, (9) 
at 

and setting X u = yj<Pj, X s - (I - Pn)X, where P N is the 
algebraic projector on X u , we have 

^+J?[ M X 1( = 0, X u (0) = P N x, (10) 

^ + M S X S = 0, X,(0) = (/ - P N )x. (11) 
at 

We set A u = [(ttpj,<Pl)} N jhsl = diagPjU^j and so, by (8), we 
may rewrite (10) in terms of y = = (x u , f*j)}j_ l as 

^-+A u y = 0, y(0) = {(P N x,<p*)} N j=i . (12) 

Since Tr(-A„) = -A\ - A 2 - ■ ■ ■ - A N < 0, it follows by Theo- 
rem 2 in [3] that there is a sequence of skew-symmetric matri- 
ces {C k }f =1 , where M = N - 1 such that the solution y to the 
Stratonovich stochastic system 

M 

dy + A u y dt = ^ C k y o d/3 k , t > 0, (13) 

has the property 

\y(t)\ < C\y(0)\e- J0 ', P-a.s., Vf > 0, (14) 

where yo > 0. The matrix C k is explicitly constructed in [3] 
and it will be used below to construct a stabilizable feedback 
controller of the form (3). Namely, we set in (3) 

N 

R k (X) = Y J C k j (x, t p*)<f >i , k=\,...,M, (15) 
where ||Cf.||f . =1 = C* 

N 

e=i 



and an are chosen in such a way that 

N 



^a i{ yej = 6ij, i,j=l,...,N. 



(17) 



Here, ytj = J tp" e Tp*d^ and since, by the unique continuation 
property (see [5], p. 157), the eigenfunction system {<p*} is li- 
nearly independent on Oq, we infer that the matrix Hy^/H^ - =J is 
not singular and, therefore, there is a unique system {an) which 
satisfies (17). Then, by (16), we see that 



(l 0o <t>i,<p*) = 6 ij , /../- 1 V. 



(18) 



As mentioned earlier, Oq is an arbitrary open subset of O. 
Theorem 1 is the main result. 

Theorem 1. The solution X to the closed loop system (4), where 
R k are defined by (15), is exponentially stable with probability 
one, that is, 

P [\X{i)\>H < Ce-y'\x\ H , W > 0] = 1, (19) 
where y > and C > are independent of co e £1 

2.1. The proof of Theorem 1 

The idea of the proof, already used in stabilization theory 
of infinite-dimensional systems with finite-dimensional unsta- 
ble subspaces, is to stabilize the finite-dimensional system (10) 
by a feedback controller of the form (3) and to reconstruct con- 
sequently the system via the infinite-dimensional stable com- 
plement (11). 

Namely, taking into account (10), (1 1), we write the closed 
loop system (4), that is, 



dX + JIXdt = P 



l 0o YjRk(X)od/3 k 



k=l 



(20) 



as 



M N 



dX u + m u X u dt = PnJ]J] C U m) ° dpk ' (21) 

k=l i,j=l 

M N 

dX s + ZKsXdt = (I ~Pn) YjTj C U ( Xu ' f*]) P((pd ° d/3k ' (22) 

k=\ !j'=l 

where X u + X s — X. 

Taking into account (18) and that X u = yj<Pj> we ma Y 
rewrite (21) as 

M 

dy e + A e y e dt = ^ C k (j yj od/3 k , 1= l,...,N, (23) 

k=l 

and so, by (14), we have that 

\ye(t)\ < Ce- y '\y e (0)\, P-a.s., Vt > 0, I = 1, ...,N. (24) 

(We have denoted by C several positive constants independent 
of t and (jj e O.) 



As regards the existence of a solution X s to (22), this is 
standard and follows from the general theory of linear infinite- 
dimensional stochastic equations. 

Now, in order to estimate X s , it is convenient to replace (22) 
by its Ito formulation (see, e.g., [10]) 



1 M 

dX s + ft s X s dt = - V [P(l 0o R k ] 2 X u dt 

L k=\ 

M 



(25) 



Taking into account that, by (7), e'^' is an exponentially stable 
semigroup on X s , without loss of generality we may assume 
that Re (!R s x, x) > y\x\^. (Otherwise, proceeding as in [6], 
we replace the scalar product (x,y) by (Qx,y) where Q is the 
solution to the Lyapunov equation J[ S Q + Q^* s = yl-) Then, 
applying Ito's formula in (25), we see that 



d\X s (f)^+2y\X,(f)^dt= (x,(f), ^[P(lo a R k ] 2 X u (t)\dt 



k=\ 
M 



(26) 



2 \P(loRk(X u ))\ 2 dt + 2 J] P(l Rk(X u ))X s d/3 k . 



k=l 



Taking into account that, by (24), |X M (0l« < e _r '|X„(0)|, we 
infer by (26) that 

E\X s {t)\lt < C\X s (Q)\ 2 e- y ' ', Vf > 0. 

Since M(t) = £ Z k= i P(lo Rk(X u )X s d/3 k is a local martingale, 
it follows by (26) and Lemma 3.1 in [6] that 

\X s (t)\-H < QX s (0)\e~ 7 ', Vf > 0, P-a.s., 
which completes the proof. 

3. The design of a real valued noise controller 

A nice feature of the feedback controller given by Theo- 
rem 1 is its simple structure. Moreover, its computation does 
not rise any special problem because the stabilizing matrices 
C k for the diagonal system (23) can be explicitly expressed. 

If all Aj, j = 1, N, are real, then the feedback controller 
{R k } given by Theorem 1 as well as the closed loop system (4) 
are real, too. However, if some Aj are complex, then the above 
feedback controller does not stabilize system (1), but its com- 
plex realization in the space H = H + iH. In this case, though 
the stabilizing feedback controller (3) is quite simple, its imple- 
mentation into real systems rises some delicate problems be- 
cause the feedback controller is in implicit form as function of 
ReX, ImX. In order to circumvent these inconvenience, we 
will derive below from the above construction a real valued 
noise feedback controller which has a stabilizing effect on the 
Oseen-Stokes system (1). 

To this end, we denote by {iffj}^ the orthogonalized system 
(via Schmidt procedure) {Re <pj, Im <pj, j = 1,...,N). Taking 



3 



into account that each <pj in the above system arises together 
with its conjugate, the dimension of the system {t//j} is again N. 
We also set 

X% = lin spanftft}* , 

and note that X™ = {Rey; y e X u ] and H = X™ ® X™, where 
Xf = {Rey; y € X s }. Moreover, the operator leaves inva- 
riant spaces X[f, Xf and if we set 

we have, of course, 

We'^'Wum ± Ce- W , Vf > 0, (27) 
for some 70 > 0. It is also easily seen that 

Tr[-^ e ] = Tr[-J?l M ] < 0. 
Then, we define as in (15) the operators 

N 

R k (x) = Y i C$j(x><Pj)*i> k=l,....,M, (28) 
ij=i 

where C k = \\C^\\fj =1 is the matrix system which stabilizes with 
probability one, via Theorem 2 in [3], the finite-dimensional 
system 

M 

dy+Al e ydt = J]c k yod/3 k . 
Here, = || if/j) \\f J=1 and ^ is given by 

N 

<t>i = ^amfig, i=l, ....,N, 
where 5% are chosen in such a way that 

N ~ c 

As in the previous case, the system {^/j} ! J =1 is still linearly inde- 
pendent on Oo which implies that det \fyij\\ + and the above 
system has a unique solution . 

Then, arrived at this point, the proof of Theorem 1 applies 
neatly to conclude that the feedback noise controller (3) with 
Rk instead of Rk is stabilizing in the probability system (9). 

Namely, we have the following stabilization result. 

Theorem 2. The solution X to the closed loop system (4) with 
Rk = Rk defined by (28) satisfies (19) in the real norm \X(t)\n- 

Remark 3. It should be emphasized that there is a close con- 
nection between the unique continuation property of eigenfunc- 
tions tpj (or <p*) of the Oseen-Stokes operator and the above 
construction of a stabilizing finite-dimensional feedback con- 
troller (3). In fact, as seen above, the design of u in the form 
(15) or (28) is essentially based on this sharp property through 
existence in the algebraic system (17). 



It should be mentioned also that, likewise the corresponding 
Ito noise controller (see [5]), the stabilizable feedback control- 
ler (3) is robust. 

Remark 4. One might speculate that the noise feedback con- 
troller 

M 

u = l 0o YjR k (X-X e )od/3 k (29) 

k=l 

inserted in the right hand side of Navier-Stokes system (2) sta- 
bilizes exponentially in probability the equilibrium solution X e 
to (2). In general, this might not be true, but it happens for Ito 
noise of the form (29) for any x in a sufficiently small neighbor- 
hood of X e (see [8]) and one might expect that the fixed point 
argument used there for the equivalent random system is still 
applicable in the present case. We expect to give details in a 
later work. 

4. Conclusions 

We have designed in this paper a Stratonovich stochastic 
feedback controller which exponentially stabilizes in probabi- 
lity a general Oseen-Stokes system from fluid dynamics. The 
controller has the support in an arbitrary open subset Oo of 
the velocity field domain O c R d , d = 2,3, and has a finite- 
dimensional linear structure which involves the dual eigenfunc- 
tions corresponding to unstable eigenvalues of the systems. The 
stabilization effect is independent of the Reynold number 1 jv 
through the dimension of the stabilizing controller (3) might 
depend on v. 
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